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simulations, we demonstrate here that the symmetry does get restored at a certain 
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allowed only if the reheating temperature after inflation is below T c . We also comment 
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1. Introduction 



It has recently been suggested that new venues for, e.g., Grand Unification, could be 
obtained by considering compact extra dimensions of an inverse radius at the TeV 
scale, with non-Abelian gauge fields propagating in the bulk (see, e.g., [TJ |2l E] and 
references therein). 

It is well known that if SU(iV) gauge fields were to couple only to adjointly charged 
matter, and the radius of a compact circular dimension is such that a weak coupling 
computation is reliable, then the system possesses a spontaneously broken global Z(iV) 
symmetry, related to the Polyakov loop in the extra direction [U El H] • This leads 
to the existence of topological domain wall defects [HI El EE HI] • 

Of course, such domain walls would carry a huge energy density, ~TeV 3 , and could 
thus not be produced under any normal circumstances. However, if the symmetry 
gets restored at a high temperature T c , domain walls could possibly be produced in a 
cosmological phase transition in the Early Universe. On the other hand, the existence 
of domain walls with energy densities exceeding about MeV 3 can be excluded experi- 
mentally, for instance, through the anisotropics in the cosmic microwave background 
radiation . This means that either particle physics models of this type are not 

realised in nature or, in analogy with how the monopole problem can be avoided in 
the case of traditional grand unified theories, the cosmological history is such that the 
reheating temperature after inflation is sufficiently below T c . 

Given the meaningful conclusions that can be drawn, it seems worthwhile to ask 
whether there really is a symmetry restoring phase transition in these models. This 
is not automatically guaranteed to be the case [T^]: there are even non-perturbatively 
established examples of broken symmetries which do not get restored at high temper- 
atures [T^]. Then one would obviously have no domain wall problem [To] 1 . 

It turns out that in the present context the issue cannot easily be solved with per- 
turbation theory, since it breaks down close to the phase transition point JH]- On the 
other hand, various non-perturbative arguments suggest that the symmetry should get 
restored JH]- Here we confirm the latter arguments with lattice simulations. 

The model we study is pure SU(2) gauge theory. We take it to live either in four 
Euclidean dimensions (4d), standing for one time direction, two space directions, and 
one compact extra dimension — or in five dimensions (5d), standing for one time 
direction, three space directions, and one compact extra dimension. A continuum limit 

1 Another possible loophole might be if the field responsible for the domain wall is very weakly 
coupled to the plasma In our case the relevant field is A y , the component of the gauge field in 
the extra direction, and it couples strongly enough even to stay in thermal equilibrium. 
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can only be approached in 4d, since a 5d pure Yang- Mills theory is not renormalisable. 
Nevertheless, for a fixed lattice spacing, the qualitative features of the phase diagram 
can be addressed also in the 5d case. 

In Sec. 121 we specify the lattice model studied, as well as the observables measured. 
We present our 4d data in Sec. El our 5d data in Sec. El and conclusions in Sec. El 



2. Basic formulation 

2.1. Lattice action and observables 

We denote the number of dimensions by d, d = 4, 5. Two of the dimensions are assumed 
compact and periodic: a time direction r = O...L T , where T = L~ l is the temperature, 
and a compact spatial dimension y = O...L y , where M = L~ l . The lengths of the 
remaining (in principle infinite) dimensions are denoted by Li, i — 1, d — 2. 

We regularise the theory by introducing a spacetime lattice, with a lattice spacing a. 
The lattice spacing is taken the same in all directions. We denote L T = aN T , L y = aN y , 
Li = aNi. The lattice action is of the standard Wilson form, 

s = /fc EEC 1 - s 1 ^)' ^ 

X fJ,<U 

where P^ u is an SU(2) plaquette, and [i, v — l,...,d. The (dimensionless) value of (3q 
determines the lattice spacing a; we return to the conversion in Sec. 13.41 

Our aim, from the point of view of the discretised theory, is to determine the phase 
diagram of the system in the three-dimensional space spanned by /3g,N t , N y . The 
spatial volume V = LiL 2 ...-Ld_2 is assumed extrapolated to infinity. 

The observables we employ in order to determine the phase diagram are the Polyakov 
loops in the r and y directions, 

I n t -i 

P T (x h y) = -Tr J] U T (x + nf), (2.2) 

* n=0 

P y ( Xl ,r) = -Tr J] U y (x + ny), (2.3) 

A n=0 

where the U T ,U y are link matrices, x = (xi,r,y), i = l,...,d— 2, and f, y are unit 
vectors in the directions of r, y. We also monitor the plaquette, (TrP^(x)). 
For both observables we measure the distribution of their volume average, 
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Let us define, in particular, 



\Pr\ 



\Pr\), X(Pr) 



VLy 
r.d-1 



((Pr) 2 )-(\Pr\)' 



(2.5) 



and in complete analogy for \P y \, x(Py)- Physically, the role of \P y \ is to tell whether 
dimensional reduction into a d — 1 dimensional effective theory takes place (\P y \ > 0) 
or not (\P y \ = 0), while the role of \P T \ is to tell whether the effective d — 1 dimensional 
theory is confining (\P T \ = 0) or not (\P T \ > 0). Let us recall that here d—1 is supposed 
to comprise the d—1 flat physical spatial dimensions, as well as the one time dimension; 
for d = 4 we are thus in effect considering a toy model of a (2+l)-dimensional world. 



2.2. Universality class 

As we vary (3c-, N T , N y , we may expect to find phase transitions. In case the transitions 
turn out to be of the second order, as is mostly the case, they are associated with some 
universality class, which we now discuss. 

We have introduced two actual order parameters for the system, P y , P T , whose phases 
are Z(2) variables. In the limit of N T ^> N y or N y ^> N T , only one of them is "critical" 
at the transition point; then the scaling goes to a good approximation according to 
the Ising model in d — 1 dimensions ^J]- For comparable N T ,N y when both can be 
critical, the universality class is that of the Z(2)xZ(2) model in d — 2 dimensions. 
The properties of this model however depend on a number of parameters, defining the 
relative strengths of the self-interactions of the two spins, as well as the interactions 
between them. In general the transition is still of the second order, but in special cases 
it can also be of the first order, at least for d — 2 = 3 20J. Below, we will indeed 
encounter first order phase transitions in the 5d case. 



2.3. Perturbative predictions 

Before turning to lattice results, let us also briefly review the predictions of continuum 
perturbation theory. The 1-loop finite temperature continuum effective potential for 
the phase of the Polyakov line, P y , has been computed in ^H], both for d = 4,5. The 
structure of the potential is such that at zero temperature, assuming weak coupling 
for the given value of M (say, M 3> A^jg in the 4d case), the Z(2) symmetry related to 
P y is broken |10j . such that \P y \ > 0. This is consistent with the fact that the theory 
undergoes dimensional reduction to d — 1 dimensions [2*Tj . 

Increasing now the temperature, the potential turns out to show no sign of sym- 
metry restoration. However, at temperatures parametrically of the order T ~ M/g 2 , 
perturbation theory is seen to break down [THj . and non-perturbative arguments can 
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N T N v Nt x N 2 



1...6 



1,2,4 



1...6 



4 



2 



1 



16 2 
16 2 

12 2 ,16 2 



Table 1: The volumes studied in the 4d case. Note that on account of the symmetry 
N T <-> N y , P T <-> P y , some of them are effectively listed twice. 

be given that the symmetry should get restored [TH] . The resulting phase diagram was 
sketched in Fig. 2 of [IHJ. 

3. Lattice results in four dimensions 

In order to study the phase diagram with lattice simulations, we shall employ various 
fixed values of N y , N y = 1,2,4. In each case, we study a number of different N T , and 
determine the phase diagram as a function of (3g- We plot the results in the plane 
(1/N T , Pg)- For a fixed (5q (fixing a), increasing 1/N T (moving right) corresponds then 
to increasing the temperature, since 1/N T = aT. 

The 4d simulation points are shown in Table d At every (3q close to a phase transi- 
tion, we study five $q- values more carefully, each with a statistics of 40K measurements 
(separated by 1 Metropolis and 6 overrelaxation sweeps). 10K sweeps are discarded for 
thermalization. Interpolation in /3a is achieved through Ferrenberg-Swendsen reweight- 
ing |221, with jackknife error estimates. The location of the susceptibility maximum is 
obtained with a parabolic interpolation through the binned, reweighted data. 

3.1. Existence of phase transitions 

In Fig. [T] we show the results obtained for \P y \, as well as the susceptibility x{Py)i as a 
function of fta, at various N T , for a fixed N y — 2. We observe that \P y \ grows rapidly 
(symmetry breaks) when /3g is increased beyond a certain point. The transition point 
is located from the maximum of x{Py)- We see that the location of the phase transition 
depends on N T . Measurements of \P T \, x(Pt) show similar phase transitions, however 
in general at some other (3q. By studying two different volumes for N y = 4 (Table [TJ), 
we know that the results are already representing the thermodynamic limit well enough 
for our purposes, as far as the position of the transition is concerned (see Fig. |21 bottom 
right). For d — 5 (Sec.EJ), we have carried out somewhat more extensive volume scaling 
tests, and they do not give rise to any concerns about our resolution. 
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Figure 1: Left: the observable \P y \ defined in Eq. ()2.5j) . for N y = 2, and various N T . 
Right: the corresponding susceptibility x(Py)i defined in Eq. (|2.5j) . For V — > oo, 
\P y \ — > for /3q to the left of the susceptibility maximum (cf. Figs. El El)- 

3.2. Phase diagrams for fixed N y 

When results such as discussed above are collected together, we obtain the 4d phase 
diagrams shown in Fig. [21 

As a first check, we may compare our results with some known limits in the literature. 
The critical values of (3c for a (2+l)d SU(2) theory at finite temperatures have been 
determined in [2B1 1211 • These results are collected in Fig. [21 top left, showing the 
transitions felt by \P T \ in the limit N y = 0, but they also determine the transition 
points felt by \P y \ at 1/N T = oo (filled boxes). The other known limit is obtained from 
(3+l)d at finite temperatures [25H2H] : these results fix the transitions felt by \P y \ at 
1/N T = (filled circles). Our data indeed interpolate between these known values. 

A second check may be obtained by considering the behaviour of the transition 
for \P y \ in the vicinity of 1/N T = 0. Indeed, the behaviour there is sensitive to the 
universality class of the system. The point 1/N T = represents a (3+l)d theory, so 
that the transition should belong to the 3d Ising universality class. As we now move 
to a small l/N T > 0, we may expect the transition point, (3 C , to change according to 
the corresponding finite-size scaling exponents. More precisely, a finite aN T acts as a 
length scale cutoff, so that if the correlation length £ critical exponent v is defined as 
£ ~ \(3 — Pd~ u , we may expect 

\(3 c (Nr) - Pc(Nr = 00)\^-^L+^ + ...). (3.1) 

1 VT ^ ' 
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Figure 2: Phase diagrams for various N y , in 4d. The lines connecting the simulation 
points are there to guide the eye only. The data for N y = at iV T > 2, as well as that 
at 1/N T = oo (filled boxes), are from the (2+l)d studies in (2S1IS|, while the limiting 
values at l/N T = (filled circles) are from the (3+l)d studies in [23 EHj For N y = 4 
we have shown results both at the spatial volume 12 2 (small crosses) and 16 2 (large 
open symbols), finding no appreciable volume dependence. 



Here we have also included u, the universal correction-to-scaling exponent. For 3d 
Ising, v ~ 0.63, uj ~ 0.8. The results, obtained by fixing the exponents u, u and fitting 
for P C (N T = oo),ci,c 2 , are shown in Fig. 01 for iVy = 2. We find that the ansatz in 
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Figure 3: Finite-size scaling of the phase transition point determined from \{Py) on 
volumes 16 2 x iV T x 2, as a function of 1/N T , together with a three-parameter fit 
according to Eq. (|3.1jl (fit for (3 C (N T = oo), ci, c 2 ). 

Eq. (|3.1|) indeed allows to fit the data well, as long as N T >N y . 

3.3. Physical implications 

After these consistency checks, we draw the following physical conclusion from Fig. |21 
Suppose we fix (3g, meaning that we fix the lattice spacing. Staying at this (3g at 
zero temperature, 1/N T = 0, let us increase M until we are in the phase where the 
symmetry related to the phase of P y is broken (region C), by decreasing N y (recall 
that M = l/(aN y )). This effectively guarantees that dimensional reduction takes 
place [2*T] . Since \P T \ = 0, the dimensionally reduced theory is in its confining phase. 
Increasing now the temperature by moving to the right (T = l/(aN T )), we observe that 
there is first a deconfinement transition to region A and then, at large enough values 
(l/iV T > 1/Ny), we do cross the transition to the symmetric phase, \P y \ = (region 
B). Thus, the phase diagram is indeed as proposed in |18j . 

3.4. Phase diagram in continuum units 

To be clearer with the comparison, we can transform the axes in Fig. |21 to physical 
units. Through dimensional transmutation, the 4d SU(2) theory is characterised by 
a physical scale which one could choose to be, say, the lightest glueball mass, the 
string tension, the Sommer scale r , or the scale parameter Aj-jg of the MS dimensional 
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regularisation scheme. Here it is convenient to set the scale instead through T , defined 
to be the critical temperature for the finite temperature phase transition in the (3+1 )d 
theory with M = 0. 

For the purposes of this paper, it is enough to carry out the conversion from /3 G , iV T , 
N y to T/T , M/T on a rather approximate level. The explicit procedure we adopt 
is as follows. For N T = 1, 16, the values of f3 G corresponding to the deconfinement 
phase transition in the (3+l)d theory have been determined in |2l)l l2*o]. Let us denote 
these values by j3 c (N T ). By a spline interpolation 2 , we treat the function (3 C {N T ) as if 
it were defined for real arguments in the interval N T = 1.0. ..16.0. We then declare that 
scaling violations (or finite a corrections) in To are small. Consequently, a given value 
of f3 G is converted to a value of aT through 

4- = N TtC = fc 1 ((3 G ). (3.2) 
Given that aN T = T -1 , aN y = M _1 , we then get 

To N T ' T N y 

Whether dependence on the lattice spacing has thus disappeared, can be checked ex- 
plicitly by carrying out the same conversions at various N y ; we do this for N y = 2,4. 
On the logarithmic scales we shall use, the residual dependence turns out to be very 
small. To address also values of (5 G larger than /3 C (16.0) = 2.74 (2H|, we run (5 G from 
there up with the perturbative 2-loop /3-function. 

The result of this procedure is shown in Fig. 0] This figure can directly be compared 
with Fig. 2 of [T%] : we find perfect qualitative agreement. It may also be noted how 
remarkably well the finite size scaling ansatz in Eq. (J3.1j) can be used to parameterise 
a fit through the data points. 

In principle it would also be interesting to check the parametric form proposed for 
the continuum critical curve (between, say, the regions A,B), T ~ M/g 2 ~ MlnM. 
We have however only a few data points, with modest resolution and without a precise 
extrapolation to the infinite volume and continuum limits, so we do not attempt that 
here. 

2 The spline interpolation introduces a free parameter, the value of the second derivative of j3 c at 
the boundary N T = 1.0, which we arbitrarily choose to vanish. On our resolution, this has little effect. 
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Figure 4: The 4d phase diagram in physical units. Filled symbols are for N y = 2, open 
ones for N y = 4. That they agree within errorbars indicates the absence of significant 
finite lattice spacing effects. The dashed curves represent the finite size scaling ansatz 
shown in Fig. El as well as its mirror with respect to the interchange r «-> y. 



4. Lattice results in five dimensions 

We then move to 5d. Considerably less can be achieved now, since the theory is not 
renormalisable and thus a continuum limit cannot be systematically approached. Nev- 
ertheless, there are a number of issues which remain the same as in 4d: the 1-loop 
finite temperature effective potential for the phase of P y is still finite and shows sym- 
metry breaking at zero temperature and small coupling ^U], but no sign of symmetry 
restoration at high temperatures jTH]. At the same time, non-perturbative arguments 
can again be presented according to which restoration should take place. We shall here 
study the issue in lattice theories defined with fixed finite lattice spacings. 

The volumes used are listed in Table El The action, observables, and simulation 
techniques follow closely those in 4d, as described in Sec. E] 



4.1. Evidence for phase transitions 

Typical signals for phase transitions in 5d are shown in Figs. El El I n contrast to 
4d, we now observe both second (Fig. EJ) and first (Fig. EJ) order phase transitions. 
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Table 2: The volumes studied in the 5d case. Note that on account of the symmetry 
N T <-> N y , P T <-> P y , some of them are effectively listed twice. The volumes in brackets 
have only been used for some of the values of N T shown. 
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Figure 5: Examples of the absolute value of the Polyakov loop, \P y \ (left), and the cor- 
responding susceptibility, x(Py) (right), as a function of fie- The absence of significant 
hysteresis on the left indicates a second order phase transition. Note that \P y \ — > for 
(3g to the left of the susceptibility maximum, as the volume is increased. The points 
on the right have been joined together by reweighting. 

The first order transition is also visible in the plaquette expectation value, as shown 
in Fig. |U1 (right). Qualitatively, the possibility of an emergence of a first order signal 
may perhaps be understood via the universality arguments presented in Sec. 12.21 

4.2. Phase diagrams for fixed N y 

The results from different N y are summarised in Fig. [7| Again the existing data from 
the (3+l)d theory [23122] allow for a first check by comparing with a number of limiting 
values (filled boxes). We find good consistency. 

A second check can be inferred from the scaling of (3 C close to 1/N T — 0, as in Sec. 13.21 
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Figure 6: Examples of the absolute value of the Polyakov loop, \P y \ (left), and the pla- 
quette expectation value, (|Tr P^) (right). The different branches show some metasta- 
bility, or hysteresis, indicating a first order phase transition. Note again that \P y \ — ► 
for Pg f° the left of the transition point, as the volume is increased. 



Now the universality class is that of the 4d Ising model, so that we expect a mean field 
exponent v = 0.5 in Eq. (|3.1|) . meaning a parabolic behaviour. In the 5d case our data 
is not precise enough to address the issue in great detail, but we do find that, say, the 
data with N T > N y at N y = 1, can be fitted reasonably well with two free parameters, 
(3 C {N T = oo) and c x (Fig.©. 

4.3. Physical implications 

Following Sec. 13.31 we can then draw the following physical conclusion from Fig. 
Suppose we fix /3 G , meaning that we fix the lattice spacing. Staying at this fie a t zero 
temperature, 1/N T = 0, let us increase M until we are in the phase where the symmetry 
related to the phase of P y is broken, such that dimensional reduction to a 4d theory 
takes place (region C). This can be achieved by decreasing N y , since M = l/(aN y ). 
Staying then on the horizontal line of a fixed /3g and increasing the temperature by 
moving to the right in 1/N T (T = l/(aN T )), we observe that at large enough values 
we may indeed cross the transition to the symmetric phase (region B). It is easy to see 
that if this happens by passing through the deconfined region A, then at the transition 
point l/N T > 1/N y , implying T c > M. Thus, the phase diagram is qualitatively as 
sketched in [18J. 

However, in contrast to the 4d case, the phenomenon described only takes place in a 



11 



N =0 N =1 




1/N 1/N 

T 1 



Figure 7: Phase diagrams for various N y , in the 5d case. The lines connecting the 
simulation points are there to guide the eye only. The values for N y — (top left), as 
well as for 1/N T = oo (filled boxes), are from j2SH2H|- Open circles and boxes denote 
second order transitions, open diamonds first order ones. 

limited range of N y , meaning that the continuum limit cannot be approached. Indeed, 
we find (Fig.0 bottom right) that for N y = 4, the phase transition close to 1/N T = is 
of the first order [2Zj, and the lines corresponding to the breaking of P T , P y have merged. 
This implies that there is no longer a low temperature phase (region C) which would 
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Figure 8: Finite-size scaling of the phase transition point determined from x{Py) on 
volumes 12 3 x N T x 1, as a function oil/N T , together with a two-parameter fit according 
to Eq. flUD (fit for f3 c (N T = oo), c h with a fixed v = 0.5, c 2 = 0.0). 

be effectively four dimensional (\P y \ > 0) and confining (\P T \ = 0) 3 . We expect this to 
remain true also for N y > 4. Since the phase diagram changes qualitatively as N y is 
increased, one cannot define a meaningful line of constant physics, where discretisation 
artifacts would gradually vanish, as correlation lengths measured in lattice units grow. 

We end with a brief comment on a related recent work, ref. [2HI- The authors study a 
(4+l)-dimensional SU(2) gauge theory, corresponding to the axis 1/N T = in our fig- 
ures, with asymmetric lattice spacings. As far as their data can be compared with ours, 
the results are consistent: with our symmetric lattice spacing, their study corresponds 
to N y <2, implying indeed a second order phase transition. Our conclusion regarding 
the existence of a continuum limit in a genuine 5d theory is however different, for the 
reasons mentioned above. Of course a continuum limit does again exist for models of 
the type in where the extra dimension is assumed inherently discretised (in other 
words, the continuum limit is only taken on the 4d planes, will the fifth dimension 
remains latticized), so that the lattice spacing is in a sense infinitely asymmetric. 

3 Since wc do not address the Higgs mechanism here, we would like the theory to be in the confine- 
ment phase even if we had in mind gauge fields related to weak interactions. 



N =1 

y 



13 



5. Conclusions 



We have discussed here the finite temperature behaviour of pure non-Abelian Yang- 
Mills theories, living in a space with a single circular extra dimension. 

If the extra dimension is small enough, dimensional reduction takes place in asymp- 
totically free (or weakly coupled) theories, as indicated by the standard perturbative 
analysis in terms of Kaluza-Klein excitations. At the same time, the circular dimension 
introduces a global Z(iV) symmetry in the system. If perturbation theory works, the 
symmetry must necessarily be broken, because the expectation value of the Polyakov 
loop in the extra direction is parametrically close to unity. 

While just introduced through perturbative arguments, the same Z(N) symmetry 
breaking can also be observed when the system is studied non-perturbatively. How- 
ever, as we have demonstrated in this paper, this symmetry gets restored at high 
temperatures, T c > M, where M is the inverse of the perimeter of the extra dimension. 
This symmetry restoration is non-trivial because it is in contrast to straightforward 
perturbation theory. 

Avoiding a domain wall problem in cosmology requires that such a symmetry restora- 
tion never takes place. We are thus lead to conclude that we would either have to change 
the model, or restrict the reheating temperature after inflation to be below T c . 

The simulations carried out to reach this conclusion were very simple. In fact, for 
d — 4, the limiting cases of zero as well as infinite radii of the extra dimension could 
be extracted from the literature. We have here simply added a number of points in 
between, utilising small scale numerical simulations. However, we believe that our 
physical conclusions, as discussed above and summarised in Fig. 01 are new. 

It is interesting to note that, at least on the resolution of a logarithmic scale, the 
phase diagram in Fig. 0] is well described by a finite size scaling ansatz in the whole 
range of interest. It is more difficult to determine the parametric behaviour of the 
continuum critical curves, but the symmetry restoration temperature is in any case 
typically above the inverse of the perimeter of the extra dimension, T c > M. 

There is less existing data available on the d = 5 case. The new points we have 
added come again from small scale simulations. On a coarse lattice, the phase diagram 
remains qualitatively the same as at d — 4. However, the continuum limit cannot be 
approached: on a finer lattice, some of the phase transitions turn into first order ones, 
and the physically relevant confining phase disappears. This means that one cannot 
define a "line of constant physics" along which to approach the continuum. Moreover, 
unlike in 4d, the correlation lengths do not diverge in lattice units as N y — > oo, if 
the first order signal continues to strengthen as expected. While the absence of a 
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continuum limit is well known, we may have elaborated on it in a somewhat new way. 

For cosmological applications, it is important to discuss the robustness of the phase 
diagram with respect to small modifications of the model. Let us reiterate that the main 
point of our study is that if a Z(iV) symmetry exists in the system and is spontaneously 
broken at zero temperature, then it does get restored at high temperatures, contrary 
to the prediction of straightforward perturbation theory. We would expect the main 
model dependence, then, to concern the zero temperature starting point, rather than 
actual finite temperature physics. 

It is certainly true that it is possible to break the Z(iV) symmetry explicitly with a 
relatively small modification of the model. While matter fields in the adjoint represen- 
tation do not break it, matter fields in the fundamental representation coupling to the 
covariant derivative D y , do break the symmetry explicitly. Another generic way may 
be to change the topology of the spacetime, for instance by orbifolding, whereby the 
zero mode of the gauge field component A y vanishes, A y = —A y = and thus P y — 1, 
or by considering more than one extra dimension of a specific structure. At the same 
time, other topological defects could of course potentially be introduced in such more 
complicated settings (see, e.g., |E1I2H])- To summarise, the existence of domain walls 
at zero temperature is by no means generic, but requires a case by case study. 

If the Z(iV) symmetry is indeed explicitly broken by some of the mechanisms men- 
tioned, then cosmological constraints are much weaker ^3] ■ Still, the existence of pos- 
sibly relatively long-lived metastable Z(iV) vacua could lead to interesting phenomena, 
in analogy with the discussion in [30J. 

Let us end by recalling that there are other cosmological constraints similar to the 
one we have discussed. In particular, the overproduction of gravitons (or more gener- 
ally very weakly coupled Kaluza-Klein excitations) at temperatures above the inverse 
radius, in case such particles propagate in the bulk, requires again a low reheating 
temperature, T<M [5T] . 
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